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Abstract 

Several quantities relevant to phenomenological studies of B^~B^ mixing are 



computed on the lattice. Our main results are fB^v Bb^ = 206(28)(7) MeV, ^ 



fBsyBBs/fsayBBa = 1.16(7). We also obtain the related quantities fBs\jBBs 
237(18) (8) MeV, /b, = 174(22)1^1^ ^eV, /s, = 204(15)™ ^eV, JeJIb, 



1.17(4)1°, fsjios = 0.74(5). After combining our results with the experimental 

world average /!\m^^^^'\ we predict Am^ = 15.8(2.1) (3.3) ps~^ . We have also com- 
puted the relevant parameters for D^-D^ mixing which may be useful in some exten- 
sions of the Standard Model. All the quantities were obtained from a quenched simu- 
lation with a non-perturbatively improved Clover action at /? = 6.2, corresponding 
to a lattice spacing = 2.7(1) GeV, on a sample of 200 gauge-field configurations. 
A discussion of the main systematic errors is also presented. 



PACS: 13.75Lb, ll.15.Ha, 12.38.Gc. 



1 Introduction 



Historically, the measurement by the UAl Collaboration of a large value for the neutral 
meson mass-difference Am^ was the first indication that top quark mass had to be very 
heavy [|I| . This mass difference is induced by B^-B^ oscillations. Since then, B^-B^ mixing 
became one of the most important ingredients of current analyses of the unitarity triangle 
and of CP violation in the Standard Model and beyond 

In the Standard Model the mixing, induced by the box diagrams with an internal top 
quark, is summarized by the next-to-leading order (NLO) formula [y, p!0[ 



q2 

= VBSoixt) \VtgVtl\'^ tubJIBb^ {q = s,d) , (1) 



where 5'o(xt) is the Inami-Lim function |jTT|, xt = m'f/M^ {rrit is the MS top mass 
mp(mp) = 165(5) GeV) and 7]b = 0.55(1) is the perturbative QCD short-distance NLO 
correction. The remaining factor, f^^iSsg, encodes the information of non-perturbative 
QCD and this is what can be computed on the lattice. 

Traditionally, the 5-parameter of the renormalized operator is defined as 

{B,\Qt''='{^^)\B,) = I mlJl BbM) ^ (2) 

where Qf^^^ = (^7^(1 — 75)9) (^7^(1 ~ 75)^') and fi is the renormalization scale. This 
definition stems from the vacuum saturation approximation (VSA) in which Bsg = 1- The 
renormalization group invariant, and scheme-independent, B parameter Bb^ of eq. (Q) is 
defined as 

Bb, = a.(/.)--«/2^" |l + BbM) > (3) 

where 70 = 4 and J depends on the scheme used for renormalizing Qf^^^ifi) (see below). 

Besides the B^-B^ amplitude, an important quantity for phenomenological applications 
is given by the ratio 

IT/ |2 

e , (4) 
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^mB^ 



where ^ = fB,\/BBjfB,\jBB,. 

In this paper, we have computed the hadronic parameters appearing in eqs. (|I|) and 
(H), using a non-perturbatively improved action, and with operators renormalized on the 



lattice with the non-perturbative method of ref. |T^, as implemented in [|13|, |1^. Our main 
results are 



/b./SI = 206(28)(7)MeV, ^= ^^^^fgl = 1.16(7) 



/srf V Bb, 



2 



2 / 



Jb^^Bb, = 237(18)(8)MeV, r,, = ^ = 1-40(18 



fB, = 174(22)™ MeV, /b. = 204(15)™ MeV , (5) 
^ = 1.17(4)+?, = 0.74(5) . 

From eq. (|T]) we write 

/\md[ps-^] = (0.153 ± 0.010) I^mI'/b.^b, [MeV^] . (6) 

We stress that what is really relevant for B^-B^ mixing, and can be directly computed 
on the lattice, is the physical amplitude, corresponding to f^^ Bb^, and not the decay 



constant /b^ and Bb^ separately. The value of fB^y Bb^ given in the first of eqs. (g) comes 
from the calculation of the amplitude, and thus it includes the correlation between the 
decay constant and the B parameter. Its final error is then smaller than that obtained by 
combining the results and errors obtained from Jb^ and Bb^ from different calculations. 
By taking \Vtd\ — ^^^J (1 — p)^ + = 0.0080(5) (where A, p, r] and A are the Wolfen- 



stein parameters) this gives 

Amrf = 0.42(12) (4) ps-\ (7) 



where the first error comes from the lattice uncertainty on fB^y Bb^ and the second from 



the error on rrit and on \ Vtd\- The most recent world average of experimental results is |jl5 

Amf^P-) = 0.473(16) ps-^ . (8) 

To predict Anig it is convenient to use the above experimental information and eq. (^ 
written as 

Am.= ^-^^e(^AmdY"' . (9) 
\Vtd\^ \mB, J 

Using our result for and \Vts\^ / \Vtd\^ ~ 1/[A^((1 - p)^ + V^] = 24.4(5.0), we get 

Am, = 15.8(2.1)(3.3)ps-^ , (10) 

to be compared with the experimental lower bound 

Am, > 12.4ps"^ . (11) 

In the above calculations, we have assumed the values of the relevant couplings, namely 
\Vtd\ and \Vts\, from the unitarity relations of the Vckm matrix. Obviously, from the 
experimental determinations of Am^ and the hadronic matrix elements we can, instead, 
constrain \Vtq\. This is what is usually done in the unitarity triangle analyses For 
the sake of comparison, we also give the values of the decay constants and their ratios from 



3 



our previous studies of these quantities, obtained using the same improved action, with a 
comparable statistics but on a larger lattice |jl6| 



fB, = 173(13)1^4MeV , /b, = 196{ll)tfMeV , (12) 
^ = 1.14(2)(1), 1^ = 0.72(2)^^^ 

JBa IDs 

The last error in the figures above represents the uncertainty coming from the extrapolation 
(linear vs quadratic) in the inverse meson mass to the B mesons. This error is not given 
in eqs. (H), because in the present study only a linear extrapolations have been made for 
reasons discussed below. 

Besides the quantities relevant to B^-B^ mixing, we also give the corresponding quan- 
tities for D^-D^ mixing. Although D^-D^ mixing is expected to be well below the expe- 



rimental limit in the Standard Model, it may be enhanced in its extensions [0. We have 
obtained 




230(14)1^ MeV , fo = 207(1 1)^^ MeV , (13) 
234(9)^1^^^ MeV , ^ = 1.13(3)1? . 



D 



Reviews of recent results of quantities considered in this paper can be found in refs. 

The remainder of this paper is as follows: in sec. |^ we give the parameters of the 
lattice simulation, illustrate the calibration of the lattice spacing and of the quark masses 
and describe the calculation of the heavy meson spectrum and decay constants; in sec. |^ 
we discuss the renormalization of the relevant operators, the calculation of their matrix 
elements and the extrapolation to the physical points; in sec. ^ we give the physical results 
and discuss the statistical and systematic errors; sec. ^ contains a comparison of our results 
with other calculations of the same quantities as well as our conclusions. 



2 Lattice Calibration and Decay Constants 

In this section we give some details about our lattice simulation and the calculation of the 
decay constants, which are a byproduct of this study. 

2.1 Lattice Setup 

The results presented in this study have been obtained on a 24^ x 48 lattice, using the 
non-perturbatively improved Clover action at /3 = 6.2 with the Clover coefficient c^^ = 
1.614, as computed in ref. |T^. The statistical sample consists of 200 independent gauge- 
field configurations. Statistical errors have been estimated by using the familiar jackknife 
procedure with 40 jacks, each obtained by decimating 5 configurations from the whole 
ensemble. The following values for the heavy- and light-quark hopping parameters, kq and 
Kg respectively, have been used: 

• 0.1352 (KgJ; 0.1349 (k^J; 0.1344 (^,3), 
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• 0.1250 (acqJ; 0.1220 (kqJ; 0.1190 (kq,) . 

Although we are not discussing hght mesons, it is important to mention some results 
which will be useful in the analysis of heavy-light meson physics, for example for the 
extrapolation/interpolation in light-quark masses. More details on the procedures used 
to calibrate the lattice spacing and fix the light-quark masses can be found in previous 
publications of our group ll6|, pO |. 



o 



The ratio of the masses of the light pseudoscalar (P) and vector (V) mesons are: 

^ = {0.597(22),,, 0.682(15),,, 0.761(8),3} . (14) 
my 



o We use the method of physical lattice planes ||21| to fix the value of the inverse lattice 
spacing Q. This method consists in the following. First we fit the vector meson mass 
to the form 

Mv{mg,mg) = ao + ai{Mp{mg,mg)f , (15) 

obtaining = 0.286(16), and ai = 1.24(13). Then we fix mv/mp to the physical 
kaon ratio {rriK* /fTT-x) and compare Mk* to the experimentally measured rriK* = 
894 MeV. In this way we obtain 

a-^ = 2.72(13) GeV . (16) 

This is the value of the inverse lattice spacing which has been used throughout this 
study. 

o For a generic physical quantity in the heavy-light meson sector, jF(mQ,mg), the 
interpolation/extrapolation in the light quark to the strange/up-down (s/d) mass is 
performed through the fit 

J^{mQ,mg) = + a'^ {Mp{mg,mg)Y , (17) 

where a^i are the fitting parameters. The light-quark mass corresponds either to the 
d-quark, when we extrapolate to Mp{md, md) = = 5.1(3) ■ 10~^ (as obtained from 
the lattice-plane method from mp/m-,^), or to the strange quark, when we extrapolate 
to Mp{ms,ms) = Mr,^^ = 0.266(17) (as inferred from eq. (|l^) by fixing Mv{ms,ms) 
to = 1020 MeV). 



2.2 Heavy-light Decay Constants 

We now discuss the heavy-light meson decay constants, which are important ingredients 
for physical predictions related to B^-B^ mixing. Since this has been extensively discussed 
in the literature (see for example ref. [jl6|), we only recall here the essential steps. 

^ As in previous publications, we use small-case letters for referring to quantities in physical units {e.g. 
mp in MeV), and capital letters for denoting the same quantities in lattice units (e.g. Alp = mpo). 
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As usual, hadron masses are extracted by fitting two-point correlation functions. For 
mesons, the standard form is 



where for J(x) = P5{x) = iQ{x)'~^^q{x), we choose t G [16,23], whereas for J(x) = 
Vi{x) = Q{x)'yiq{x), we choose t G [19,23]. The time intervals are established in 
the standard way [i.e. after inspection of the corresponding effective masses). The 
resulting masses of the pseudoscalar and vector heavy-light mesons, with the light 
quark interpolated/extrapolated to s/d (see eq. (p^), are listed in tab. |l]. We also 
checked that on the same interval (t G [16,23]), the masses of pseudoscalar mesons 
extracted from the correlation function C^p(t) {Ao{x) = Q{x)jo"y5<l{x)), are indeed 
the same. 

In physical units, the masses directly accessed from our simulation, are: 
mp, = {1.75(8) GeV, 2.02(9) GeV, 2.26(11) GeV} , 
mp^ = {1.85(7) GeV, 2.11(9) GeV, 2.38(10) GeV} . (19) 

The pseudoscalar meson (P) decay constant is defined as 

{0\Ao\Pip = 0)) =iFpMp , (20) 

where the 'hat' symbol means that the appropriate renormalization constant Za is 
taken into account. To determine the decay constant, one computes the ratio 

EdMimm) M^j,,„hfM,(I-t)) . (21) 



When working with improved Wilson fermions, the axial current is improved by 
adding to the operator d^P with a suitable coefficient 

(0|io|P) = Za({0\Ao\P) + CA{0\adoP5\P)) = iMp{fP + CAaPj.'^) , (22) 



where the constant, ca = —0.037, was computed in ref. [^. It is actually easy to 
see that 

.(1) _ v^i 



= • (23) 

The size of the correcting term CAaFp^ / Fp^ for our data is in the range (5 -r- 7)%. 

As far as the normalization constant is concerned, its improved value is given by 
(choice 1.) 

ZAiniQ, m,) = Za{0) (1 + ^^am) , (24) 
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where m = [rriQ + m^)/2. In the above equation, we have taken 



(25) 



The value of the critical hopping parameter, Kc = 0.13580(1), is obtained from the 
linear fit Mp ^ rriq — 0, whereas the value of improvement coefficient 6^ = 1-24 in 
( P^D is taken from boosted perturbation theory |^ Q For a recent non-perturbative 
determination of Ba, see ref. . . 



In our numerical estimates we have also used 



ZA{mQ,mg) = Za{0) 



1 



am 



[1 + Ba am) 



(26) 



This equation (choice 2.) differs from the choice 1. in that it contains higher order 
tree-level mass corrections through the so-called KLM factor [^. As for the value of 
Za = Za{0) in the chiral limit, we take the nonperturbative determination Za{0) = 
0.80 from ref [^. This value has been also non-perturbatively computed in refs. 
{Za = 0.82) and {Za{0) = 0.81), using different approaches. The values of the 
lattice decay constants Fp for different heavy-quark masses are listed in the third 
column of tab. ^ 



"Flavor" content {Q — q) 


Mp 




Mv 


Q3- s 


0.866(4) 


0.0884(21) 


0.898(4) 




0.832(5) 


0.0781(35) 


0.869(7) 


Q2- s 


0.777(4) 


0.0858(20) 


0.812(4) 


Q2-d 


0.742(3) 


0.0761(29) 


0.782(6) 


Qi- s 


0.681(5) 


0.0853(19) 


0.721(5) 


Qi-d 


0.643(3) 


0.0755(24) 


0.690(5) 



Table 1: Mass spectrum of heavy-light pseudoscalar and vector mesons, with the hght- 
quark masses extrapolated/interpolated to the d/s quark, using eq. (p!?]). The meson 
decay constants are also given. All quantities are in lattice units. 



To estimate the values of the physical quantities related to S-mesons {ms^ = 5.28 GeV, 
and ms^ = 5.38 GeV), the extrapolations to the physical heavy-quark masses necessarily 

^ By working with three light quark species only, we were unable to make a quadratic fit in the quark 
masses as done in ref. [|o| . 
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207(ll)t;j tl MeV 


234(9)l|j tl MeV 


J Us 1 J 

1.13(3)t!] 


174(22)t^ 1^ MeV 


204(15)1^ tl MeV 


1.17(4)1° t; 



^0 

0.46(7) 


-0.73(13) GeV 


0.56(5) GeV^/^ 


$('^=^) 
-0.81(7) GeV 



Table 2: Summary of physical results for the decay constants. The parameters of eq. 
are also listed. Our value for $i agrees well with QCD sum rule estimates, namely $i 
-0.9(2) GeV m. 



rely on the heavy-quark symmetry. Accordingly, the decay constants scale as 

/p=5^!^fn-^^ + ...). (27) 

Apart from an overall (mild) logarithmic correction, the coefficients $(mp) and $'(mp), 
are non-perturbative quantities which we extract from a fit to the lattice data. Since we 
work with three heavy quarks (three values of kq), we have only included the leading 1/mp 
correction by fitting the decay constant to the expression 



= *« (^1 + ^ (28) 

In tab. ^, we give the physical results, and the values of parameters $o,i; iii physical 
units. The statistical errors are given in parentheses. The systematic errors require some 
explanation: the ffist one is obtained from the differences between results obtained by 
using Za with choice 1. or 2; the second is obtained by comparing the central values of 
the procedure described above to those extracted from the study of the ratio Fp/My. 

We have not included in the systematic errors: i) the uncertainty due to quenching. 
Pioneering attempts for estimating quenching errors [^-[0 show an increase of about 
20 ^ 30 MeV of the value of the decay constants, ii) the error due to the truncation of the 
1/mp expansion in eq. (pSj) . The latter is usually estimated by considering also a quadratic 
term (cx 1/mp). With only three heavy-quark masses, we have not done such a fit. From 
our previous experience |T^, we know that this error has a negligible effect on the value 
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of the D-meson decay constants, whereas for the B meson it represents a major source of 
systematics (~ 30 MeV). 

Finally, we present an estimate of J'b^ obtained using the experimental measurement of 
= 241(32) MeV 0. With the ratio fs/ fos obtained in our simulation, table |], we get 




(latt.) 

= 0.74(4)+^ +2 = 0.74 ±0.05 



/b = 184 ± 24(exp.) ± 12(theo.) MeV (29) 



3 Matrix Elements of the Renormalized Operators 

In this section, we discuss the renormalization of the relevant four-fermion operators 
and describe the calculation of the matrix elements introduced in eq. (H). Since we are 
interested to the B^-B^ mixing amplitude, only the parity-even part of the operator, 
Q = {h^q)(b-f^q) + (bl^nm)(hi^'l^Q)^ will be discussed in the following. 

3.1 Renormalization of the Four-fermion Operators 

On the lattice, the renormalized operator Q{^) takes the form 

Q{lx) = Z{ii,gl) [o, + j2^i{gl)0^ , (30) 

where the Oi denotes a bare operator and we choose to work in the following parity-even 
basis 

03 = bqbq -b-f^qb'j^q (31) 

04 = bqbq + }r)^qb-i^q 

05 = ba^^q ba^^q . 

To obtain the renormalized operator, two steps are then necessary. First, one has to 
subtract operators of the same dimension, which (on the lattice) mix with Oi (this is the 
consequence of the presence of the Wilson term in the fermion action, which explicitly 
breaks chiral symmetry). This means that one has to compute the subtraction constants 
Ai{gQ). After an appropriate subtraction of the operators Oi, with i ^ 1, the matrix 
elements still need (to be finite) an overall renormalization constant provided by Z{^,gQ). 
A detailed discussion of the mixing matrix can be found in ref. [|T^. In this work, we use 



the method of refs. [|T3|, |T^ to calculate the matching and the renormalization constants 
non-perturbatively. The computation, at three different values of the scale {fia = 0.71, 
1.00, 1.42), is performed in the RI-MOM renormalization scheme, in Landau gauge. The 
results are listed in table ^ For sufficiently large /i, we can use the available NLO anomalous 
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Scale a/i 




A2 


A3 


A4 


As 


0.71 


0.630(11) 


-0.052(3) 


-0.016(2) 


0.014(2) 


0.001(1) 


1.01 


0.583(7) 


-0.055(2) 


-0.017(2) 


0.016(2) 


0.002(2) 


1.42 


0.601(8) 


-0.069(2) 


-0.022(1) 


0.016(1) 


0.007(2) 



Table 3: Summary of the results for renormalization and subtraction constants (eq. 
evaluated non-perturbatively at (3 = 6.2 at three different scales /i. Z[fi,gQ) is the renor- 
malization constant in the RI-MOM scheme. 



dimension of the operator {Q{fi))^^ mom [33|^ define the renormalization group invariant 
(RGI) operator: 

{Q)= w-\^^){Q{^^)) ^ a.(/i)-^°/2^° |l + {Q{^^)) . (32) 

where 

17397 - 2070ri/ + lOAn} ^ ^ 

70 = 4 ; Jri_mom = 6(33 - 2n,y + ^ ^ " ^^^^ 

We also give the NLO part relevant in the MS scheme 

13095 - 1626/2/ + 8n} 
^ 6(33 - 2n/)2 • ^^^^ 

By using ^^(m^) = 0.118, and setting rij = 4, we obtain 

/i = {1.9 GeV, 2.73 GeV, 3.85 GeV} , 

^(/^)m-mom = { 1.408, 1.453, 1.489 }. (35) 

In quenched calculations, there is always the embarrassing question whether to use the 
physical value of ag, with the complete formulae for the anomalous dimensions and the 
/3-function, or a "quenched value" of the coupling constant (which suffers from intrinsic 
ambiguities), together with anomalous dimensions and the /3-function computed for n/ = 0. 
The present case is particularly lucky, since by computing w(/i) in the unquenched (n^ = 4) 
and quenched cases (with Aq^cd = 250 MeV), one finds values of Bb^ which differ by less 
than 2%. 

We give in passing, w{mh)^ = 1.477, which will be used later on for a comparison 
with the results of other authors. In the perturbative calculation of w{mh)^, we have 
used TTift = 4.25 GeV 
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3.2 Matrix elements 



We have computed the relevant three-point correlation functions: 

C§\p,q;tp,,tp,) = I dxdye'^P^y-^^-^^H0\P,ix,tp,)Qi6,0)^^^ . (36) 

When the sources are sufficiently separated and far away from the origin, the lowest pseu- 
doscalar mesons are isolated and one has 



C§)(p,g;tp„tpJ {PimiPiP)) ^^^-'"''''''^^ (37) 

where Zp = {0\P^\P). In order to eliminate irrelevant factors and to extract the B para- 
meter we form the following ratio 

where /i indicates the scale dependence. The two point correlation functions are both 
multiplied by the axial-current renormalization constant in the chiral limit 2^a(0). We have 
not used the improved renormalization constant ZA{fnQ,mq) since we have not attempted 
to improve the operator Q. This is beyond the scope of the present study and would 
demand the inclusion of many operators of dimension seven, whose coefficients are presently 
unknown. Our hope is that a part of the uncertainties of 0{rha) cancel in the ratio (^). 
This has been computed at |p| = |g| = 0, corresponding to 

) ^ 3(m^M,B,(,,. (39) 

In this way we obtain the main results of this work. We also considered reasonably low 
momentum injections to the sources but the additional noise makes these data not useful 
in practice. For this reason these data will be ignored in the following. 

In order to show the quality of the signal, in fig. we show the ratio TZ{fi) defined in 
eq. ( pHD at /i = 3.8 GeV {fia = 1.42), for a specific combination of the heavy and light 
hopping parameters (kq = kq^ = 0.1220 and Kq = Kg^ = 0.1349). The ratio is plotted 
as a function of the time-distance of one of the two sources. The other has been fixed to 
tp^ = 16. We also performed the simulation by fixing tp^ = 12 and tp^ = 20. From the 
study of the two-point correlation functions, one isolates the lightest pseudoscalar state 
around t = 16. That makes tp^ = 12 rather small for our purposes. For tp-^ = 20, the 
sources are not far enough from each other because of the periodic boundary conditions, 
and this may spoil the signal. After inspection of the ratios TZ{fi) (for each combination 
of the heavy and light quarks), a stability plateau is found in the range 28 < t < 33. From 
this plateau we extract the values of TZin), that is the B parameters, for each combination 
of quark masses. At this point, we extrapolate the light quark to the s- and (i-quark mass 
by fitting the B parameter to eq. (|1^. This is illustrated in fig. ^ for kq.^ = 0.1220 at 
/i = 3.8 GeV. In tab. ^, we give a detailed list of results for the RI-MOM Bp^{fi), as well 
as for the RGI Bp^ (see eq. (|3^) ). 
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1.00 



^ 0.80 



0.60 



|I = 2.7 GeV 



+ + + + + + + 



20 



30 

Time (t) 



40 



Figure 1: Ratio 7^(yu) at fia ^ 1.42, wiiicii corresponds to fi c::^ 3.8 GeV. Tie values of the 
heavy and hght hopping parameters are Kg = 0.1349 and kq = 0.1220, respectively. 



0.90 



1^ 



0.85 



«3 



0.80 



0.05 



0.1 



Figure 2: Fit of = Bp^{fj,) in the light-quark mass according to eq.(|^). Empty 

symbols denote the data obtained from our simulation ; filled symbols represent Bp^ln) 
and Bp^{fi), obtained after interpolation to the strange and extrapolation to the down 
quark mass. The heavy quark corresponds to kq^ = 0.1220. 
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Scale (/i) 


1.9 GeV 


2.7 GeV 


3.8 GeV 


Light quark 


q = s 


q = d 


q = s 


q = d 


q = s 


q = d 


Bp^{fi;KQ = 0.1250) 
Bp^{kq = 0.1250) 


0.853(16) 
1.201(23) 


0.842(23) 
1.186(32) 


0.789(15) 
1.148(21) 


0.779(21) 
1.132(30) 


0.831(15) 
1.237(21) 


0.821(21) 
1.222(29) 


5p,(/i;KQ = 0.1220) 
5p^(kq = 0.1220) 


0.875(19) 
1.233(26) 


0.867(27) 
1.222(39) 


0.810(18) 
1.178(25) 


0.803(25) 
1.167(37) 


0.851(18) 
1.264(25) 


0.846(25) 
1.259(38) 


5p,(/i; Kg = 0.1190) 
5p,(«:o = 0.1190) 


0.880(18) 
1.239(25) 


0.879(29) 
1.238(41) 


0.815(17) 
1.184(23) 


0.814(27) 
1.182(39) 


0.855(16) 
1.273(24) 


0.854(27) 
1.272(39) 



Table 4: i?p^ for the three values of kq and with the hght-quark mass extrapo- 
lated/interpolated to the d/s quarks using eq. (|17|). The results are given for both the 
RI-MOM scheme and the RGI case. 



It is interesting to check whether the scale dependence of the extracted Bp^{fi) is well 
described by the NLO anomalous dimension (eq.( |32|)). In fig. we plot the evolution 
of i?pj~^^^(/i). The three curves (dashed, full and dotted) have been normalized to the 
B parameter at the scale fia = 0.71, 1.01 and 1.42 respectively, as given in table ^ We 
note that the point at yua = 1.01 is sensibly lower than what expected on the basis of the 
perturbative evolution. This effect is a consequence of a fluctuation of the value of Z{fi) at 
this value of the scale as shown in table |^. We checked that the fluctuation is enhanced by 
the extrapolation of Z{n) to the chiral limit and consider it as part of the uncertainty in the 
determination of the matrix element. The central values of Bp^ quoted below correspond 
to the result obtained converting the RI-MOM B parameters at fia = 1.42 to the RGI one. 
The difference with the results obtained by using the other two scales, namely fia = 1 and 
/ia = 0.7, is included in the estimate of the systematic uncertainty. 
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Figure 3: NLO evolution of Bf-^^^{fi) (kq = 0.1220;. The dashed, sohd and dotted 
hnes are obtained by starting the evolution from the lattice result (filled circles) at fia = 
0.71, 1.01, 1.42, respectively. 



4 Physical results 

In this section, we derive the physical amphtudes from the B parameters computed at 
fixed niQ. Before discussing the extrapolation in the heavy-quark mass, let us summarize 
our results. We have computed the matrix elements of the operator Q{fi) in the RI-MOM 
scheme at three different scales fi. The dependence on light quark is rather smooth, and 
we have easily performed the extrapolation/interpolation to the d and s quarks. By using 
the NLO evolution, we have checked the stability of our results. For each value of the 
heavy-quark mass, from the B parameter in the RI-MOM scheme we have computed Bp^ 
by using the result at ~ 3.8 GeV. This is the only step where perturbation theory comes 
in the game. At this point we dispose of the following values: 

kq = { 0.1250 , 0.1220 , 0.1190 } 
Bp^ = {1.2A{2nZ l-26(2)tn^, 1.27(2)^^^ 

Bp^ = {1.22(3)tn^, l-mnZ l-27(4)t°:°° } (40) 

In order to extrapolate the above results ito the b quark, we rely on the heavy-quark 
symmetry as in the case of decay constants and fit the 13 parameters to the expression 

Bp. ^c-(l.|i). (41, 

The result is shown in fig. ^, and the main numbers are listed in tab. |. We stress again 
that we are not able to include terms of order 1/mp in our fit, because we work with 
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three values of the heavy-quark mass only. 



Note, 



however, that our results for the slope 
55], where the heavy quark was treated 
non-relativistically. We will comment on this in the next section. These results, when 



Ci' are much smaller than the predictions of ref 



















1.24(4)1°:°° 


1.38(11)1°:°° 


1.35(5)_o:o8 


0.98(5) 



Co 

1.46(2) 


11=^) 

Cl 

-0.28(16) GeV 


Co 

1.41(6) 


Cl 

-0.24(6) GeV 



Table 5: Physical results for 13p^ parameters obtained in this work. In addition, we give 
the parameters Cq^i from eq. (0). 

combined with those in tab. |^, give the following final numbers 



fB,\/BB, = 206(28)(7) MeV, /b.V^b. = 237(18) (8) MeV ; 



JbA B 



Bp, , , n / mo 



2 \ (cxp.) 



1.16(7) , r,d = e\^\ = 1.40(18) . (42) 



m 



In practice, for the reasons discussed in the introduction, we have extracted the value of 



fs^y Bb^-, as well as the other quantities appearing in the above equation, directly from 
the calculation of the mixing amplitude. 

For completeness, and for comparison with other calculations, we also give the B pa- 
rameter in the MS scheme 

5f(m,) = 0.93(8)1°:°, 

B^f{m,) = 0.92(3)1°:°. (43) 
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Figure 4: Filled circles correspond to the B-meson sector. We see the 1/mp negative slope 
for and Bb^ ■ Third picture shows the SU(3) breaking for B-parameter, which is clearly 
compatible with 1.0 for all heavy-light mesons. 



5 Comparison with Other Calculations and Conclu- 



sions 



In this section we compare our results with those obtained with different lattice calculations 
or with other methods. 



On the lattice, two approaches have been used so far to compute B^-W mixing: the 
direct one (as done in this work) which extrapolates in the heavy-quark mass and the one 
where the heavy quark is treated with an effective theory. In the latter, either the Heavy 
Quark Effective Theory (HQET) or Non-Relativistic QCD (NRQCD) were used. There is 
a large number of lattice studies in which the HQET was invoked. Unfortunately, most 
of them were plagued by a mistake in the calculation of the renormalization constant. 
Only recently, this has been corrected in ref . |^ , and the data of refs. P7| , have been 
reanalyzed 0. On the other hand, very recently, the NRQCD treatment of the heavy 
quark has been employed to compute Bb^ |^|- The authors also calculated the relevant 
0{l/mp) corrections. The slope in 1/mp that they observe for B^^{mh) is roughly a factor 
of 3 larger than ours (or the one by ref. ||3^). We stress that there is a general argument 
which shows that the uncertainty on the 1/mp corrections can be easily as large as the 
corrections themselves unless perturbation theory on the leading term is not pushed to very 



high orders This argument found confirmation by the explicit calculation of the 0{al) 



perturbative corrections to the heavy-quark mass in the HQET [|ll. Moreover, even the 



lowest order result of ref. |^ in 1/mp is suspicious, since they combine the renormalization 
constants computed in the HQET with the computation of the matrix elements of NRQCD. 



^This is the reason why we quote results from these three papers only. 
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Figure 5: Comparison of our results with those of ref. as a function of the inverse 
heavy-meson mass in lattice units. 



Until the renormahzation of the relevant operators will not be completed in a consistent 
theoretical framework, the results of ref. ||3^ should not be used for comparisons with other 
approaches. 

In the direct approaches, a full relativistic treatment is given to both heavy and light 



quarks |18|. In this case a major source of systematic error comes from the extrapolation 



of the results to the physical B mesons. Recent results for B^-B^ mixing can be found 



in refs. and In the second paper, an extrapolation to the continuum (a 0) 

obtained using non-improved Wilson fermions, has been attempted. They obtain a large 
central value for the ratio r^^ (r^^ = 1.76(10)142), albeit with large errors. The central 
value is difficult to reconcile with the value of ^ and of the decay constants produced 
by the same authors with the same set of data. In the present work and the one by 
UKQCD ||39|, improved fermions are used but without extrapolation to the continuum. 
The differences between the present study and that of ||3^ is that in our case the action is 
improved non-perturbatively whereas UKQCD uses the mean-field improved action, and 
that we renormalize non-perturbatively the relevant four-fermion operators. The overall 
agreement is excellent. This is true also at the values of the heavy-quark masses were 
lattice measurements have been actually made as shown in fig. ^ Their extrapolated value 
has an error smaller than ours, probably because they also use a point corresponding to a 
small value of mp, for which the heavy-quark expansion may be questionable. In table 
for comparison, we list several recent results of the direct and the HQET approaches. 
For completeness, let us mention that also QCD sum rules were employed to compute 
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Reference 




BbJBb^ 


This work 


i.38(ii)t:°° 


0.98(5) 


m 


1.41(6)^:°^ 


0.98(3) 


m 


1.41(18) 


~1 


36| + |37| 


1.29(8)(6) 




36| + |38| 


1.26(6)(6) 





Table 6: Recent lattice results on Bp^ (at NLO accuracy). The upper part of the table 
refer to results obtained by extrapolation to the B-mesons, whereas in the lower part, the 
HQET was used. The SU(3) breaking ratios are also listed (when available). 



Bb^ @- Radiative corrections were also included in ref. where the value for B^^{mb) 
was found to be compatible with 1 (VSA value). The most recent QCD sum rule esti- 
mate has been given by Chernyak who quotes a value lower than ours and closer to 
the results obtained from the HQET at lowest order in l/mp: 5^^^(5.04 GeV) = 0.82, 

corresponding to Bb^ = 1.26. The effect of SU(3) breaking in BbJBb^, was studied in 
the framework of the HQET, combined with chiral perturbation theory ||44]. Their result 
suggests that SU(3) breaking for the Bb parameter is practically negligible as confirmed 
by many lattice calculations. The approach has been extended to quenched chiral pertur- 
bation theory [Q, where a pessimistic estimate of quenching errors (~ 10-20% for the 
matrix element) was quoted. 
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